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1.  INTRODUCTION 


Experimental  results  are  cumulative  if  \>i  aggregate  they  unify  and  extend  empirical  relations  and 
theoretical  structures  that  may  be  obscur.a  in  individual  investigations.  Empirical  cumulativeness,  which 
Hedges  (1987)  describes  as  "...  the  degree  of  agreement  among  replicated  experiments  or  the  degree  to 
whidi  related  experimental  results  fit  into  a  simple  pauem  that  makes  conceptual  sense,"  is  the  focus  of 
this  report.  Glass  (1976)  was  among  the  first  to  recommend  the  use  of  quantiutive  procedures  in 
integrative  research  reviews  and  to  introduce  the  term  "meta-analysis"  to  cover  the  collection  of  such 
IHOcedures.  Meu-analysis  claims  certain  classical  statistical  procedures,  as  well  as  ^roaches  develt^ied 
specifically  for  research  synthesis,  and  has  found  at^lication  in  the  social  and  bicdogical  sciences.  The 
unification  of  experimental  results  obtained  by  differem  investigators  (grating  independently  with  their 
own  experimental  protocol  and  sometimes  using  different  methods  of  analysis,  is  the  kernel  of  meta¬ 
analysis.  A  comprehensive  treatment  of  this  subject  is  given  by  Hedges  and  Olkin  (1985). 

Measurement  in  the  physical  sciences  is  generally  regarded  as  highly  accurate,  and  although  some 
variability  is  inevitable,  the  variation  itself  is  thought  to  be  insignificant  from  a  practical  standpoint. 
Counterexamples  to  this  notion  are  {dentifiil.  even  in  carefully  conducted  experiments.  Consider,  for 
instance,  the  situation  described  by  Touloukian  (1975)  involving  two  sets  of  measurements  taken  on  the 
thermal  conductivity  of  gadolinium.  These  data,  shown  in  Figure  1, "...  are  for  the  same  sample,  measured 
in  the  same  laboratory  two  years  apart  in  1967  and  1969.  The  accuracy  of  curve  1  was  stated  as  within 
1%  and  that  of  curve  2  as  0.5%  ..."  and  yet,  the  curves  differ  by  more  than  several  hundred  percent  at 
higher  values  of  temperature. 

Physical  scientists  normally  bring  a  careful  qualitative  analysis  to  their  research  studies.  If  prudently 
employed,  interrogative  statistics,  which  are  part  of  meta-analysis,  have  a  contribution  to  make  in  the 
physical  sciences  as  well. 

After  data  have  been  collected  according  to  a  carefully  constructed  experimental  design  (e.g.,  see 
Montgomery  1991)  the  main  reason  for  determining  a  correlation  (a  regression  analysis)  is  to  examine  the 
effects  that  some  variables  exert,  or  appear  to  exert,  on  others.  Even  when  no  intuitive  physical 
relationship  is  ^rparent,  regression  analysis  may  provide  a  convenient  summary  of  the  data.  The  summary 
can  be  accomplished  in  a  number  of  ways  and  has  been  an  active  area  of  investigation  since  the  time  of 
A.  M.  Legendre  (1752-1833),  who  published  the  first  account  of  regression  by  least  squares  in  1805. 
Section  2  of  this  report  reviews  the  correlations  that  have  been  advanced  for  steady  flow  through  inert 
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Figure  1.  Thermal  conductivity  of  gadolinium. 

spherically  packed  beds  and  some  of  the  consequences  of  the  atteitdant  data  analysis.  In  Section  3,  a 
meta-analysis  of  the  gas  flow  resistance  measurements  is  undertaken.  Section  4  contains  a  summary  and 
main  conclusions. 

2.  REGRESSION  ANALYSIS  OF  GAS  FLOW  RESISTANCE  MEASUREMENTS 

Ergun  (1952),  Kuo  and  Nydegger  (1978),  and  Jones  and  Krier  (1983)  have  proposed  models  relating 
coefficient  of  drag  to  Reynolds  number  for  steady  flow  through  packed  beds  of  inert  s|^res.  However, 
the  correlations  were  developed  under  different  experimental  regimens.  Robbins  and  Gough  (1978)  also 
investigated  coefficient  of  drag  at  a  high  Reynolds  number  but  presented  their  results  in  terms  of  a  frictitm 
factor  fj'  =  F,4*^c/(1  -  4>)].  which  is  the  ratio  of  coefficient  of  drag  F^,  and  Reynolds  number  Re  scaled 
by  a  solids  loading  factor  (1  -  41). 

In  comparing  Ergun ’s  relation 


F^  =  150  +  1.75 


Re 


1  -4» 


(1) 


2 


to  that  of  Kuo  and  Nydegger 


Fv  =  276.23  +  5.05 


f  Re 


^0.87 

> 


(2) 


or  of  Jones  and  Krier 


Fv 


150  +  3.89 


f  Re 


(3) 


a  slight  notational  difference  portends  substantial  complications.  Equation  1  is  a  simple  linear  model. 
Equations  2  and  3  are  nonlinear  in  the  sense  that  one  or  more  parameters  appear  nonlinearly.  Nonlinearity 
compl:  ates  the  statistical  arudysis  of  the  data  since  determining  qjpropriate  dioices  for  the  parameters 
in  Equations  2  and  3  becomes  a  computationally  intensive  optimization  procedure,  and  inference  about 
the  resultant  relation  and  parameters  becomes  much  more  tentative.  The  mathematical  underpinnings  of 
nonlinear  regression  will  not  support  as  much  in  the  way  of  statistical  inferertce  or  hypothesis  testing  as 
is  available  for  linear  regression,  hi  general,  nonlinear  models  should  be  avoided  unless  there  is  a 
compelling  reasrai  for  their  use.  Draper  and  Smith  (1981)  discuss  this  issue  in  greater  detail. 

Standard  regression  procedures  are  developed  under  several  assumptions.  Fundamental  among  these 
is  that  the  le^nse  (here,  Fy)  is  measured  with  error  but  the  predictorfs)  (here.  Re  and  are  measured 
without  error.  Jones  and  Krier  provide  estimates  of  error  for  Fy,  Re,  and  confirming  that  this 
assumption  is  not  met,  and  call  into  question  the  efficacy  of  the  resultant  correlations.  Sometimes  an 
attempt  to  circumvent  this  requirement  is  undertaken  by  arguing  that  the  error  in  predictor  measurement 
is  sufficiently  small  as  to  be  ignored  when  compared  to  the  range  of  the  predictor  variable.  If  this  claim 
is  invoked,  reliance  upon  any  resultant  representation  must  be  tempered  accordingly. 


Since  a  correlation  provides  a  convenient  represenution  of  the  available  data,  a  direct  attempt  at 
evaluating  the  adequacy  of  a  regression  equation  involves  an  examination  of  the  differences  between  the 

measurements  taken  and  the  values  predicted  by  the  equation.  These  differences,  Fy.  -  Fy.',  i  =  1, 2 . n, 

are  called  residuals;  Fy.  is  an  experimentally  determined  value  of  drag  coefficient,  and  Fy.'  is  the 
corresponding  value  predicted  by  the  regression  equation.  A  residual  plot  for  Equation  3  is  shown  in 


Figure  2.  These  plots  may  serve  as  a  diagnostic  tool  in  addition  to  assessing  the  adequacy  of  a  fitted 
regression  model. 

Figure  2  strongly  suggests  that  another  crucial  regression  assumption  is  not  satisfied.  The  variance 
of  the  residuals  does  not  ^pear  constant  over  the  range  of  Re'  =  Re/(1  -  ^);  moreover,  the  departure  from 
the  fitted  equation  is  systematic  with  bead  diameter.  D|,.  Jones  and  Krier  recommend  reverting  to  the 
relation  (Equation  2)  proposed  by  Kuo  and  Nydegger  to  describe  their  own  measurements  taken  for  6-mm 
beads.  This  recommendation  is  data  specific  and  is  difficult  to  justify  in  general.  They  conjecture  that 
an  interaction  between  bead  size  and  tube  diameter  may  be  present,  but  this  requires  quantitative 
substantiation.  In  general,  weighted  least  squares,  or  a  transformation  of  the  observations  F^.  before 
regression,  are  potential  corrective  procedures  suggested  by  this  residual  pattern. 


Figure  2.  Residuals  vs.  particle  diameter  D^:  Jones  and  Krier  data  with  6-mm  beads  excluded. 

Nonlinear  regression  algorithms  normally  seek  to  minimize  the  sum  of  the  squared  residuals — as  in 
ordinary  linear  regression — ^in  attempting  to  determine  the  "best”  choice  of  parameters  to  model  the  data. 
These  procedures  have  previously  been  cited  as  computationally  intensive.  More  specifically,  they  are 
iterative  and  may  diverge  or  converge  to  local  extrema,  depending  upon  the  choice  of  initial  conditions. 
Through  a  systematic  selection  of  initial  conditions,  tlie  authors  determined  that  the  equation 
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F^  =  61+2.7  .  (4) 

provides  an  improved  representation  of  the  data  reported  by  Jones  and  Krier. 

The  root  mean  square  error  (RMSE).  an  estimate  of  the  standard  deviation  of  the  residuals  and  a 
commonly  used  measure  for  adequacy  of  fit.  is  reduced  by  20%  ctxnpared  to  that  corresponding  to 
Equatitm  3.  The  measurements  taken  on  the  6-mm  beads,  the  chief  contributor  to  heterogeneity  of 
variance,  have  been  excluded  horn  the  regression,  making  the  comparison  with  Jones  and  Krier  direct 
A  reduction  of  ore-fifth  in  RMSE  is  not  by  itself  a  stunning  imivovement  but  it  does  focus  more  sharply 
on  the  underlying  physical  process.  The  residual  plot  for  Equation  4  still  exhibits  the  uivlesirable  pattern 
of  undercover)  fitting  categories  of  bead  diameter  but  is  an  improvement  compared  to  the  display  in 
Figure  2. 

The  data  collected  by  Robbins  and  Gough  (1978.  1979).  which  "...  correspond  to  several  tests 
performed  on  several  occasions"  for  beds  of  q)heres.  right  circular  cylinders,  and  multiperforated  cylinders, 
may  be  transformed  into  units  appropriate  for  comparison  through  the  relationship  f/  &  F^4^e/(1  -  ^)]. 
The  authors  confined  the  analysis  to  data  taken  on  1.27-mm-diameter  lead  shot  and  on  4.76-mm  and 
7.94  mm-diameter  steel  stores,  and  determined  the  equation 

F^  = -237 3.14  .  (5) 

U  ~^) 

for  representation  of  flow  through  sidierically  packed  beds.  Equations  4  and  5  are  shown,  along  with  the 
previously  established  correlations  (1  and  2)  in  Hgure  3. 

Transforming  the  variables  (Re'.  F^)  by  taking  logarithms,  which  was  suggested  by  the  residual  plot 
in  Figure  2.  effectively  linearizes  the  data.  In  regression  analysis,  a  measure  of  precision  of  the  regression 
line  which  is  used  in  addition  to  RMSE,  is  given  by  a  statistic  denoted  as  R^.  R^  assumes  values  in  the 
unit  interval  [0, 1]  and  quantifies  the  amount  of  variation  in  the  response  accounted  for  by  the  regression 
line.  Values  close  to  one  are  highly  desirable,  indicating  that  the  regression  has  effectively  accounted  for 
most  of  the  variation  in  the  re^nse.  The  regression  line  determined  after  logarithmic  transformation  of 


5 


Figure  3.  Proposed  models  for  relating  coefficient  of  drag  and  Reynolds  number. 

the  Jones  and  Krier  data  has  =  0.98.  Tte  transfoimed  Robbins  and  Gough  data  have  R^  =  0.99.  These 
values  are  so  close  to  1.0  that  pursuit  of  a  nonlinear  model  is  difficult  to  justify  mathematically. 

Comparison  between  linear  models  and  nonlinear  models  is  difficult.  RMSE  values  cannot  be 
compared  across  the  transformation,  and  a  well-defined  R^  statistic  for  nonlinear  models  does  not  exist. 

3.  META-ANALYSIS  OF  GAS  FLOW  RESISTANCE  MEASUREMENTS 

Consider  in  aggregate  the  correlations  that  have  been  advanced  for  gas  flow  resistance  measurements 
through  si^rically  packed  beds.  For  the  nonlinear  models,  a  statistical  resampling  plan  is  applied,  whose 
goal  is  to  extract  information  from  a  set  of  data  through  repeated  inspection.  The  procedure  is  called  the 
"bootstrap,"  named  to  convey  its  self-help  attributes,  and  it  tries  to  address  an  important  problem  in  data 
T«iuction — after  an  estimate  of  some  parameter  is  computed,  what  accuracy  can  be  attached  to  the 
estimate?  Accuracy  here  refers  to  the  "±  something”  that  often  accompanies  statistical  estimates,  and  may 
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be  conveyed  through  such  devices  as  variance.  RMSE.  or  confidence  interval.  For  the  log-linear  model, 
the  available  data  are  directly  combined. 


The  authors  are  hindered  in  fully  exploiting  a  meta-arudysis  ai^xoach  by  the  inability  to  obtain  all  the 
pertinent  experimental  data.  It  is  unfortunate  that  experimental  data  ate  tx)t  routinely  archived  after 
collection;  otherwise,  additional  information  that  it  may  hold  is  lost  to  extraction  by  subsequem 
investigations  and  by  alternate  statistical  methods.  The  data  of  Jones  and  Krier  and  of  Robbins  and 
Gough,  were  accessible.  With  these  data,  this  repon  proceeds  as  far  as  statistical  prudence  permits. 

3.1  Bootstrarwwng  Regression  Correlations.  Detailed  descriptions  of  the  bootstrap  and  accounts  of 
its  successful  af^lications  are  amply  documented  (e.g.,  Efron  [1979, 1982],  Efron  and  Tibshirani  [1985], 
LePage  and  Billard  [1992]).  The  computational  contrivartce  that  the  bootstrap  procedure  exploits  is  the 
generation  of  perturbed  data  sets  from  a  single  set  of  data  through  sampling  with  replacement  Specific 
to  this  study,  the  set  of  paired  observations  taken  on  coefficient  of  drag  and  Reynolds  number,  ((Fy^. 

RejO . (Fy^.  Re„')).  that  is  the  basis  for  a  reported  correlation,  is  sampled  with  replacement  to  generate 

another  set  ((Fy  *,  Re/*) . (Fy  *,  Re/*)}  whose  elements  arc  copies  (with  duplication)  of  the  original 

measurements.  This  set  is  called  a  bootstra(^)ed  data  set  The  process  of  sampling  with  lei^acemau  to 
generate  bootstrapped  data  sets  is  repeated  many  times. 

If  a  correlation  is  determined  for  each  bootstrapped  data  set  and  its  equation  plotted,  a  sense  of  the 
sensitivity  of  the  regression  line  to  perturbation  of  the  original  data  comes  into  focus.  In  Figure  4,  the 
results  of  1,(X)0  replications  of  this  process  ate  pictured.  The  outermost  lines  indicate  boundaries  within 
which  the  correlation  (5)  might  be  expected  to  lie  if  the  original  data  set  were  simply  perturbed.  They 
were  obtam  I  from  the  maxima  and  minima  of  the  drag  coefficient  predicted  for  particular  values  of 
Re'.^  The  envelope  constructed  for  correlation  (5)  contains  correlation  (4).  This  suggests  that  no 
significant  difference  between  these  empirical  relations  exists.  Similar  results  are  otxained  if  we  begin 
with  correlation  (4);  correlation  (5)  will  lie  within  die  corresponding  confidence  envelope.  Ctxisideration 
of  perturbed  data  is  highly  tqipropriate  here,  since  experimental  results  cannot  be  expected  to  be 
reproduced,  even  if  the  experiment  is  replicated  under  tightly  controlled  conditions.  The  theoretical 
justification  for  the  use  of  bootstrapped  data  is  given  by  Efron  (1982). 


^  More  precisely,  the  values  represent  extreme  quantiles  after  all  the  FyS  have  been  ranked;  their  values  are  not  essentially 
different  from  maxima  and  minima. 
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Figured  Bootstrapped  confidence  envelopes  for  nonlinear  regression  (based  oti  Robbins-Gourii 
data). 

The  relationship  of  Kuo  and  Nydegger,  for  which  the  experimental  data  was  not  accessible,  was 
deteimined  for  a  single  diameter  bead,  =  0.83. 

3.2  Log-linear  Regression.  Figure  5  displays  the  logarithmic  transfoimed  data  of  Jones  and  Krier, 
and  Robbins  and  Gough,  combined.  The  fitted  line  for  these  data  is 

logF^  =  0.49  +  0.90  logRe'  ;  (6) 

included  in  the  regression  are  the  data  taken  on  6-mm  beads,  which  were  previously  excluded. 

Visually,  the  data  appear  linear  after  transfonnatitm.  Statistically,  the  R^-value  for  the  regression  is 
0.99,  making  the  fitted  line  a  highly  satisfactory  representatitm  of  these  data  for  iL  practical  purposes. 
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combined 


4.  SUMMARY  AND  CONCLUSIONS 


-  Jones  and  Krier.  Robbins  and  CJoueh  data 


For  Reynolds  numbers  exceeding  10^,  a  more  effective  representation  and  data  analysis  than  presently 
availatde  can  be  obtained  after  logarithmic  transfonnation  of  the  data.  This  linearizes  the  data  and 
removes  the  necessity  for  nonlinear  regression  techniques.  The  equation 

logF,  =  0.49  +  0.90  logRe'  (7) 

is  an  effective  description  of  the  available  experimental  data. 


If  a  represoitation  of  the  form 


Re 
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is  required,  then  Jones  and  Krier’s  results  are  more  effectively  reflected  through  tire  equation 

F,=  61  +  2.7|y^J  .  (9) 

and  Robbins  and  Gough’s  data  restricted  to  ^rhericaily  packed  beds  provide  the  relation 

F^  = -237 +  3.14  .  (10) 

^1 

but  here  again,  approximate  confidence  envelopes  ctxistructed  with  the  aid  of  the  bootstrap  suggest  that 
these  relations  can  be  combined  without  loss  of  underlying  i^iysical  insight.  In  total,  the  statistical 
analysis  sui^rts  the  combination  of  the  various  correlations,  for  the  stated  test  conditions,  into  a  single 
relationship. 

While  it  is  quite  reasonable  to  suspect  an  interaction  between  the  geometry  of  tube  and  packing, 
perhaps  reflected  through  the  ratio  more  extensive  testing  is  required  to  establish  this  relation. 
Hopefully,  this  will  be  done  in  accordance  with  a  formal  statistical  experimental  design  to  minimize  testing 
and  maximize  extraction  of  information. 

G.  E.  P.  Box,  an  important  contemporary  statistician,  has  remarked  that  "No  model  is  correct,  but 
some  are  usefttl."  In  this  spirit  these  remarks  are  offered,  along  vidth  the  hope  for  an  iiKremental  move 
toward  a  more  useful  model. 
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NOMENCLATURE 


Db 

Dc 

V 


Re' 


u 


Pi.  i  =  0.  1.  2 


AP 

P 

4> 


(1  -♦) 


P 


spherical  particle  (bead)  diameter 
test  chamber  diameter 
Fy[Re/(l  -  ♦)],  friction  factor 


pu 


^2 


1  -  ^ 


,  coefficient  of  drag 


i-th  observed  value  of  the  drag  coefficiem 

predicted  drag  coefficient  corresponding  to  the  i-th  observed  value 

length  scale 

Rep  ^  =  puDi,^/p.  Reynolds  number 
Re/(1  -  ♦) 


Reynolds  number  based  particle  size 

average  gas  velocity 

model  coefficient 

change  in  pressure 

density 

porosity  of  the  packed  bed 
solids  loading 
gas  viscosity 
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